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Abstract Small and elongated, cool and dense blob-like structures are being reported with 
high resolution telescopes in physically different regions throughout the solar atmosphere. 
Their detection and the understanding of their formation, morphology and thermodynami- 
0\ ' cal characteristics can provide important information on their hosting environment, espe- 

cially concerning the magnetic field, whose understanding constitutes a major problem in 
solar physics. An example of such blobs is coronal rain, a phenomenon of thermal non- 
equilibrium observed in active region loops, which consists of cool and dense chromospheric 
blobs falling along loop-like paths from coronal heights. So far, only off-limb coronal rain 
(~| ' has been observed and few reports on the phenomenon exist. In the present work, several 

C^h! datasets of on-disk Ha observations with the CRisp Imaging SpectroPolarimeter (CRISP) 

at the Swedish 1-in Solar Telescope (SST) are analyzed. A special family of on-disk blobs 
\^ I is selected for each dataset and a statistical analysis is carried out on their dynamics, mor- 

"^ . phology and temperatures. All characteristics present distributions which are very similar to 

C^ ' reported coronal rain statistics. We discuss possible interpretations considering other similar 

blob-like structures reported so far and show that a coronal rain interpretation is the most 
likely one. Their chromospheric nature and the projection effects (which eliminate all direct 
^ ] possibility of height estimation) on one side, and their small sizes, fast dynamics, and espe- 

f~^ ■ cially, their faint character (offering low contrast with the background intensity) on the other 

C^^ I side, are found as the main causes for the absence until now of the detection of this on-disk 

^^ . coronal rain counterpart. 

rn 

pg , 1 Introduction 

The determination of the coronal magnetic field constitutes a cornerstone in solar physics. 
The response of the plasma, its structure, dynamics and thermodynamic properties largely 
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depend on the structure and strength of the coronal magnetic field. The difficulty in assessing 
the latter is thus a major hindrance in t he general understanding of the pla s ma ch aracteristics 
in the solar atmosphere. As shown in I Antolin and Rouppe van der VoortI (l2012h . one of the 
most attractive features of coronal rain is that, due to the very small sizes involved, it can 
act as a probe of the local magnetic field structure and strength, as well as give valuable 
information about the local thermodynamic conditions inside loops. 

Coronal rain is a phenomenon of active region coronae. It corresponds to cool and dense 
matter and not waves JDe Groof, Berghmans. and van Driel-Gesztelvil2004l ; lDe Groof et al.l 
l2005h that are observed oflf-li mb falling from co r onal heights down to the lower solar atmo- 
sphere along loop-like paths (lKawaguchilll97d : iLerovll 19721) . Typically observed in chro- 
mospheric lines such as Ha or Ca ii H, similar features, but with absorption p rofiles in EUV 
spectral lines, have been observed frequently in vyarm active region loops ( Foukall 197a. 



J978':'Schriiver'200l': loi^Shea. Baneriee. and Dovlel2007l:lTOpathi et al.l2069l : IUgarte-Urra. Warren, and Brooksl 

2009; Kamio et al. 2 01 ih . If these reports do correspond to the same phenomenon, it sug- 

gests a sc e nario in which coronal rain is rather common in active region loops. On one hand, 

ISchriiven KOOlh . using various channels of TRACE, has estimated the occurrence rate of 

coronal rain in active region loops to be at most once every two days, suggesting a sporadic 

character for the phenomenon. On the othe r hand, in Ha observa tions at the limb with the 

CRisp Imaging SpectroPolarimeter (CR ISP, Scharmer et al.l2008|), at the Swedish 1 -m Solar 

Telescope (SST. IScharmer et al.ll2003ah . lAntolin and Rouppe van der VoortI ( 120121) present 

a ubiquitous character of coronal rain, and show that it is composed of a myriad of small 

blobs, with sizes that are, on average, 300 km in width and 700 km in length. Furthermore, 

if close enough together and in large enough quantities, the blobs are seen as large clumps 

termed 'showers'. Such large events occur sporadically and can have widths up to a few 

Mm, thus suggesting that what has been observed with coarser resolution instruments such 

as TRACE or SOHO/EIT actually corresponds to these kinds of showers. 

Observational papers on coronal rain show a broad distribution of falling speeds, with 
values up to 120 km s"' or more, and average values around 60-70 km s"U. The respec- 
tive downward accelerations are in general lower than those due to the effective grav- 
ity along loops, suggesting the presence of other forces, for instance from the increas- 
ing gas pressure at lower atm ospheric heights. In Ca ii H observations with Hinode/SOT, 
lAntolin and Verwichtd ( l201ll) detected transverse MHD waves through coronal rain track- 
ing and showed that the pressure force from the waves can also explain the low falling 
speeds. By measuring the wave characteristics, they further estimated the magnetic field 
strength through coronal seismology techniques. 

The temperatures that have been attributed so far to the phenomenon range from tran- 
sition region down to chromospheric temperatur es, according to the lines in which it has 
been observed, from Ne vii 465 A down to Ha JLevine and Withbroe 1977; Miiller et al 



12005 ). By using the width of the line in Ha as a proxy, lAntolin and Rouppe van der Voort 



1 2OI2I) give upper limits for Ha coronal rain centered around 7000 K but with a long tail 



to higher temperatures up to 5 x lO'' K. It has not been possible to measure densities di- 
rectly so far, but estimates are made based on numerical simulations of the phenomenon 
JMuller. H ansteen, and Petej2003l;lMuller, Peter, an d Hansteen |2004l;lTsiklauri et al.l2004 

Imc 



iMendoza-Briceno Sigalotti. and Erdelvil2005l ; lMok et al.i2008l ; 



Antolin. Shibata, and VissersI 



' Most of these observations are done with imaging instruments and thus the speeds conespond to 
projected values in the plane of the sky. However, since all reports correspond to off-limb observa- 
tions, the values are close to the t otal v elocities, as confirmed by the spectropolarimetric observations in 
lAntolin and Rouppe van der VoortH2012h . 
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l20ld : lMurawski. Zaqarashvili. and Nakariakovll201lh . These produce bl obs with densitie s 
at 10'" - lO" cm"^, values that are similar to those found in prominences (IHiravamall 1985r) . 

Coronal rain is part of a general phenomenon of thermal instability in plasmas that 
takes place when ever radiation losses locally overcome the heating input in a specific struc- 
ture ( lFieldlll965h . This can happen following a density perturbation (anything leading to 
an increase in the density, such as a shock wave) since the cooling increases faster than 
linearly with density. As a consequence, temperature and pressure drop in the perturbed 
region, accreting gas from the surroundings and forming an increasingly larger conde nsa- 
tion, which can also be seen as an entropy mode (Murawski, Zaaarashvili. and Nakari akovl 
l201lh . This cascading eflFect proceeds until heating and cooling balance again at some 
lower temperatures and higher densities (and when pressure balance is regained). First 
applied to prominences (|Parkenll953c lKleczeklll957h . the thermal instability phenomenon 
has also been suggested for being responsible of observed structure at larger scale s: plan- 
etary nebulae ( IZanstralll955l) . spiral arms condensing out from the galactic halo JSpitzen 
Il956h . co ndensation of interstellar clouds (lFieldlll962h . filamentary structure in interstellar 
medium llCoxlll972r) . Numerical simulations in the last 30 years have largely contributed 
to the understanding of this phenomenon (Goldsmith 1971 i Hildnei 1 1974 ; iMok et al.lll990l: 



Antiochos and Klimchuklll99l] : iDahlburg. Antiochos. and Klimchuklll99a ; lAntiochos et al.l 



1999l : lMok et al.ll2008l) 



In the case of the solar corona, this thermal instability is known as 'thermal non-equilibrium' 
or also 'catastrophic cooling'. The high densities necessary for the instability onset are 
thought to be achieved through footpoint heating. The latter consists of having the heat- 
ing concentrated towards the footpoints of coronal loops. In such a scenario, chromospheric 
evaporation together with direct mass injection into the corona from the heating events en- 
sure a dense corona. As thermal conduction is insufficient in transporting enough energy to 
the dense corona, the coronal temperature is consequently reduced over time. This is en- 
sured by the negative slope of the radiative loss function in the coronal temperature range. 
At some point a critical state is reached in which any density perturbation (for instance, any 
shock wave traveling along the loop) is enough to trigger the fast regime of this thermal 
instability. The catastrophic cooling and condensation that ensues implies recombination of 
elements. The partially ionized clumps form and become visible in cool lines. Depending 
on the existing forces (gravity, magnetic and gas pressure gradients), the clumps either fall 
(coronal rain), or remain suspended (prominences). Since the material can remain in the 
corona over long periods of time supported by the magnetic field, the element population, as 
well as the thermodynamics, can differ significantly between coronal rain and prominences. 
It is thus important to distinguish both phenomena. 

The importance of coronal rain is further stressed through its link to coronal heating, be- 
ing the observational signature of thermal non-equilibrium (and hence to footpoint heating). 
As shown by the numerical simulations cited previously, thermal non-equilibrium is very 
sensitive to parameters such as the heating scale height and the loop length. By studying 
corona l rain we can then learn about the heating mechanisms. lAntolin. Shibata. and VissersI 
(120101) showed that Alfven wave heating, a strong coronal heating candidate, is not a pre- 
dominant heating mechanism in loops with coronal rain. When propagating from the pho- 
tosphere into the corona, Alfven waves can nonlinearly convert to longitudinal modes due 
to density fluctuations, wave-to-wave interaction, and deformation of the wave shape dur- 
ing propagation (Vasheghani Farahani et al. 201jJ). These modes subsequently steepen into 
shock s and heat the plasma uniformly along the loop ( IMorivasu et al.l2004ljAntolin and Shibatal 
|2010|) . thus avoiding the loss of thermal equilibrium in the corona. 
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As Stated above, recent high spatial and temporal resolution instruments such as SST/CRISP, 
Hinode/SOT, and SDO/AIA are presenting a scenario in which coronal rain is a rather com- 
mon phenomenon of active regions. However, all coronal rain studies so far have been per- 
formed off-limb, where emission profiles allow an easy detection of the blobs, less hindrance 
from projection effects exists and thus heights above the surface and trajectories are more 
easily determined. A natural question is then wh ether coronal rain can still be observed fur - 
ther onto the disk. A first attempt is presented in I Antolin and Rouppe van der VoortI l l2012h . 
where, apart from the ofF-limb cases, on-disk blobs with absorption profiles in the same 
dataset are also analyzed and shown to be part of the same phenomenon through dynam- 
ics, sizes and temperature comparison. It is the aim of this work to extend the study in the 
same direction and, through Ha observations with SST/CRISP at various on-disk locations 
(different datasets), present the on-disk counterpart of coronal rain. 

The work is organized as follows. In Section|3]the various datasets are presented, as well 
as a description of the data manipulation and reduction. In Section[3]we explain the methods 
used for the determination of the blobs characteristics. Results are presented in Section |4] 
followed by discussion in Section[5]and finally conclusions in Section|6l 



2 Observations and reduction 

2.1 Observational setup 

In this paper we consider five data sets obtained in the Ha line with the CRISP instrument at 
the SST located at La Palma (Spain). The CRISP instrument is a dual Fabry-Perot interfer- 
ometer (FPI) capable of tuning to the desired line position within a spectral line in < 50 ms. 
As the light comes down from the telescope tower, it first passes through an optical chopper 
and a wavelength selection prefilter, where the latter is mounted on a filterwheel that can 
change to the required prefilter position within 1-2 s. After the prefilter, a few percents of 
the light is branched off to a camera that serves as a wide-band reference channel for the 
image post-processing, while most of the light is guided through the FPI and a polarizing 
beam splitter before falling onto two additional cameras. The exposures obtained with these 
high-speed low-noise Samov CAMIMIOO CCD cameras, which all three expose at a rate 
of 35 frames per second with an exposure time of 17 ms, are synchronized by means of the 
optical chopper. 

In combination with the high spatial resolution of the SST, CRISP thus allows for ob- 
servations at high spatial, spectral and temporal resolution, which is necessary to effectively 
observe the small and fast scale dynamics of the solar atmosphere. Although for a major- 
ity of data sets considered in this study multiple lines were observed, in the following we 
only consider the Ho- data (for which CRISP has a transmission FWHM of 6.6 pm and its 
corresponding prefilter has a FWHM of 4.9 A). 



2.2 Data acquisition 

The first data set was obtained on 10 June 2008 between 7:26-8:00 UT, at a 4.2 s cadence. 
The Ha line was sampled at 10 equidistant line positions ranging from -0.4 A to H-0.5 A 
around line centre at 0.1 A spacing and 4.2 s cadence. The field-of-view (FOV) covers an 
area of 67x67 arcsec", showing AR 10998 (classified the day after as such) and contains 



On-disk coronal rain 



a small growing sunspot and surrounding plage located roughly at {x,y) ~ (-847,-115), 
corresponding to // = 0.43. 

Data set 2 was obtained on 11 June 2008 between 7:55-8:32 UT at a 6.2 s cadence and 
samples the Ho- line at 21 equidistant positions out to ± 1.1 A. The data cover roughly 
the same FOV (both in area size and in target) as data set 1. The small sunspot displays a 
rudimentary penumbra surrounded by bright patches of network and is located at (jc, y) x 
(-684, -152), corresponding to ^=0 .67. This data set has be en studied before, both for Eller- 
man b ombs in lWatanabe et al.l ( 120 llh and flocculent flows in lvissers and Rouppe van der VoortI 
ll2012h . 

The third data set was obtained on 27 June 2010 between 13:31-13:58 UT as part of 
a dual-line observing program where nearly co-temporal line scans were obtained in both 
Ha and Ca n 8542 A, resulting in an overall cadence of 17.1s. The Ho- line scans consist of 
41 equidistant line positions spaced at 85 mA and sampling the spectral line out to ±1.7 A 
around line centre. The FOV area covers 56x56 arcsec^ and contains a large sunspot with a 
well-formed penumbra at a viewing angle of //=0.25. 

In the fourth data set, obtained on 28 June 2010, the same active region as on 27 June 
(classified as AR 11084 by then) was observed between 08:16-09:08 UT, but now as part 
of a triple-line observation program including (in addition to the same lines as observed the 
day before) also full Stokes measurementsat -48 mA in the Fe i 6302 A line. The Ha line 
scans were extended with respect to 27 June to cover 45 line positions and thereby covering 
the range out to ±1.9 A around line centre, still at the same equidistant spacing of 85 mA. 
The FOV, though containing the same sunspot (this time at a viewing angle //=0.64 on the 
East side of the Sun), is slightly smaller in area size (54x53 arcsec-), while with 22.4 s the 
cadence for this data set is slightly longer than on 27 June. 

The fifth and last data set, obtained on 6 July 2010 between 08:02-09: 10 UT, covers 
again the same active region AR 11084, at a viewing angle of //=0.61 on the West side of 
the Sun. The observing program run on that day was the same as on 28 June: A triple-line 
program, where the Ha line is sampled at 45 equidistant line positions at an overall cadence 
of 22.4 s, only covering an even smaller FOV area of 52x52 arcsec^. 

For data set 1 and 2, the camera pixel size is 0.071 arcsec px"', while that for data sets 3 
through 5 is 0.0592 arcsec px"' . In both cases this falls well below the Rayleigh diffraction 
limit for Ha at the SST, which lies at 0.17 arcsec. Table [T] gives a summarized overview 



Table 1 Overview of the Ha CRISP data sets analysed in this study. 



Data set 


Target 


Date 


Location 


A 


AA 


At 


Duration 








(,") 


[A] 


[mA] 


[s] 


[min] 


1 


AR 10998 


10 Jun. 2008 


0.43 


-0.4-0.5 


100 


4.2 


34 


2 


AR 10998 


1 1 Jun. 2008 


0.67 


-1.1-1.1 


100 


6.2 


37 


3 


AR 11084 


27 Jun. 2010 


0.25 


-1.7-17 


85 


17.1 


27 


4 


AR 11084 


28 Jun. 2010 


0.64 


-1.9-1.9 


85 


22.4 


52 


5 


AR 11084 


06 Jun. 2010 


0.61 


-1.9-1.9 


85 


22.4 


68 



of the data set properties described above, while Figure [T] shows samples of the observed 
regions in Ha. 
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(a)SST/CRISP 10 June 2008 07:37:12 UT 
Ha line center Ha + 0.4 A 




60 


» 






55- 




^^^^H^Hi 


i 


" ^° 






K 


a 






^L>' 


X 45- 






P 


40- 




"f y 


< 


35- 


1 


.^m^^ - 





40 45 50 55 

X [arcsecl 



(b) SST/CRISP 11 June 2008 08:19:33 UT 

Ha line center Ha + 0.7 A 
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(c) SST/CRISP 27 June 2010 13:39:17 UT 
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(d) SST/CRISP 28 June 2010 08:37:40 UT 
Ha line center Ha + 0.9 A 
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(e) SST/CRISP 6 July 2010 08:58:28 UT 

Ha line center Ha - 0.8 A 
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Figure 1 Sample images from the data sets analyzed in this study showing good quality stills at moments 
with coronal rain occurrences. Full FOV Yla Hne center images are displayed in the left-hand panels, while in 
the right-hand panels show zoomed-in portions of the FOV (indicated by the white box in the left-hand panels) 
at a selected offsets from Une center (specified at the top of each panel) where some blobs are best observed. 
In all panels (parts of) the paths tracing the coronal rain blobs (see also Section im have been overlain. The 
white arrows in the right-hand panels indicate the location of some of these blobs at that particular time. 



2.3 Reduction 

Although multiple lines were observed in three out of five data sets, only the Ha data has 
been considered in this paper. We shall therefor only focus on the reduction of those data in 
the descriptions given in this Section. The quality of the data obtained with CRISP at the SST 
is already greatly improved during the data collection proc ess through the use of real-time 
tip-tilt correction and the adaptive optics system at the SST ( IScharmer et al.l2003bl) . In addi- 
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SST/CRISP Ha +0.8 A 10-06-2008 



SST/CRISP Ha+ 0.8 A 10-06-2008 
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Figure 2 Images illustrating the methods used for calculating the length (left) and width (right) of a blob 
along its trajectory, belonging to the 10 June 2008 dataset. The left panel is a space-time diagram where 
the spatial coordinate indicates the distance along the blob path, and the time coordinate indicates the time 
from the beginning of the observation. The right panel is an image of the same blob (indicated by the black 
arrow in that image) at a specific time. The white dots in both figures (contours on the right panel) show 
the extension of the respectively calculated blob size (length and width) along the fall. Only the successful 
size measurements are shown. Piece-wise segments are fitted along the trajectory thus defining measurements 
(^i,ti)b (see Section|3]for details). The extension of one such measurement is indicated by the white arrows 
in the left panel and by the black contours in the right panel. In the overlaying plots we show in dotted lines 
the intensity of the background at a specific time ( (indicated by the black an'ow in the left panel), r^'^ ° 
(in arbitrary units in the plot) where ,v denotes distance (in Mm in the plot) along the blob's path or across it, 
depending on whether it is a length or width measurement. In solid lines we show the intensity of the blob at 
that time with subtraction of the background, I^'^ ' (its position along the path is shown by the black arrow 
in the right panel), and in dashed fines the result of the gaussian fit to the latter The FWHM of the gaussian 
fits are taken as a measure of the size (length and width) of the blob. For the corresponding measurement we 
have a length of 1 . 17 ± 0.04 Mm and a width of 0.34 ± 0.02 Mm. 



tion, M ulti-Object Multi-Frame Blind Deconvolution (MOMFBD. Ivan Noort. Rouppe van der Voort. and Lofdahll 
l20Q5h , further improves the quality of the observations in the image post-processing stage 
by removing most of the remaining high-order seeing effects. For that procedure, each im- 
age (i.e., at each line position within the scan) is divided into 64x64 px- overlapping sub- 
fields that are processed as a single MOMFBD restoration. Precise alignment of the re- 
stored narrow-band CRISP images is achieved using the wide-band anchor exposures. Fur- 
ther details concerning the MOMFBD post -processing of similar data sets can be found in 
Ivan Noort and Rouppe Vtin der Voord (12008!) . Additional data post-processing includes cor- 
rection for the wide-band prefilter transmission profile, image rotation as a r esult of the 
alt-azi muth mount of the SST, as well as destretching of the images following IShine et alj 
( Il994l ). in order to remove most remaining small-scale rubber-sheet seeing effects. 
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3 Methods 

For the sa ke of comparison with coronal rain s tudies , especially with the recent observational 
study by lAntolin and Rouppe van der Voorj ( l2012h . in this work we analyze the dynam- 
ics (velocities and accelerations), sizes (lengths and widths) and thermodynamic properties 
(temperatures) of the observed blobs in an attem pt to assess their nature. For this purpose w e 
employ the same methods as in the analysis by lAntolin and Rouppe van der Voord ( l2012h . 
but provide more detailed explanations here. 



3.1 Dynamics 

Tracing of the blobs is performed with t he help of the CRisp SPectral EXplorer (CRISPEX, 
Ivissers and Rouppe van der Voora2012h and its auxiliary program TANAT (Timeslice ANAl- 
ysis Tool), both widget-based tools programmed in the Interactive Data Language (IDL), and 
which enable the easy browsing and analysis of the image and spectral data, the determi- 
nation of blob paths, extraction and further analysis of space-time diagrams. Fitting piece- 
wise segments along the blobs trajectories in the space-time diagrams we can estimate the 
projected velocities (in the plane of the sky). With this method, the standard deviation for 
projected velocities is estimated to be roughly 5 km s"' . 

Along a trajectory s defined by a blob, several other blobs can be observed falling down. 
Let Bs be the set of all blobs observed for a trajectory s. The fitting of piece wise segments 
along a trajectory naturally defines a family of measurements or lines (.«,-, r,)^, / = 1, ..., w^ 
for each blob b e B^, where nt is the total number of measurements for that blob along the 
trajectory. The procedure described in the previous paragraph assigns to each measurement 
/ = 1, ..., Ill, a projected velocity v^^.^.. We take /','q = (I^q ) as the observed spectral profile 
of the blob corresponding to the i'-th measurement, where the angle notation denotes the 
mean over the points defined by the line (xi, ti)t,. In order to calculate the Doppler velocities 
of a blob along its trajectory the spectral profile of the background needs to be subtracted 
from the observed spectral profile of the blob. This is first achieved by determining the set of 
times Ts at which the trajectory s does not present any blob-like structure. The background 
spectrum corresponding to each measurement (.«,-, f,)^ (/ = 1, ..., w^) is then defined as Z^'^^, = 

{I^Q ''), such that f; + to e t, and where the angle notation denotes the average over the 
line (xi, ti + to)/,. That is, we calculate the mean intensity over equivalent lines in the space- 
time diagram as those corresponding to the measurements for the blob, but at times without 
blobs. In this way, since the timescale in which the background changes is normally larger 
than the timescale in which a blob falls, we ensure a more accurate background spectrum. 
The true spectral profile of the blob is then defined as /^'j = I'^.^^ - /^'q. Since the blobs 

are observed on-disk we have absorption spectral profiles, and we ensure /^' j > at the 
wavelength position where the blob is best observed. In order to reduce errors, only clearly 
discernible blobs are selected, and only measurements for which I'^'^ is above a specific 
threshold. The Doppler velocity of a blob along its trajectory v'^^ j^^, assigned to the i'-th 
m easurement is calculated usirig the fi rst moment with respect to wavelength, a method used 
bv lRouppe van der Voort et al.l (120091) to calculate the velocities of the disk-counterparts of 
type II spicules: 
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where c is the velocity of light, Aq the wavelength at line center, and where the integration 
range is set by the minimum and maximum wavelengths for which I'^ ^ > 0. In order to 
estimate the error involved in this calculation we performed the integration varying both end 
points by various amounts ell: /imm + dl, Amax - dl. The Doppler velocity with this method 
was then set as the mean over the resulting values, and the standard deviation gives us an 
estimate of the error. 

The obtained velocities were checked with two other methods. The first involves a single 
gaussian fit of /',' j . The wavelength interval where the gaussian fit is made is basically the 
same as the [/Imin, /(max] range defined above and is allowed to vary in the same way in 
order to estimate the errors involved. The Doppler velocity with this method is taken as 
^Do ier2 ~ pC'^gauss — '^o), where /(gauss is the maximum of the gaussian fit. The second 
check is simply v'^^ j^^, = j-(/tmax - •^o), where /(max is the location of the maximum of 
/^' J in the range [/Imin, AjaaA- While all three methods deliver similar results we note that the 
gaussian method generally delivers slightly lower velocity magnitudes, and the maximum 
method normally delivers larger values. 

Having the projected and Doppler velocities we can then calculate the total velocities 

corresponding to the i'-th measurement: vj^^ = Jiv'l^jg^^ + (^Donnier^"- ^^^ resulting accel- 
eration for the i'-th measurement is then simply defined as a'-* = Av'^^^IAt'^, where the 
diff^erence is calculated between two consecutive measurements and where the t\,^ is the 
time at the midpoint of a measurement {xi, ti)b. Given that the projected velocities are in 
most cases larger than the measured Doppler velocities, the error in the acceleration mea- 
surement is mostly dependent on the projected velocity measurement. We have estimated 
the error on the acceleration measurement to be roughly 0.03 km s~^. 



3.2 Sizes 

For calculating the length /'•* of a blob at the /-th measurement we perform a gaussian 
fit to the intensity profile along the blob's length at the wavelength position of maximum 
blob intensity and at a particular time. This is illustrated in the left panel of Figure |3] More 
precisely, given the blob position {x, t) G (jc,, f,);, we first determine the wavelength /(,■ at 
which the blob intensity is largest. We then select a specific length interval /, covering the 
blob's length at time t. The background spectrum is removed from the observed spectral 
profile, as explained in the previous Section, obtaining the true profile of the blob for a 
specific wavelength /(, and time t: I'^ j = I'^ ^^ - 1'^ q, which only depends on position. Next, 
we determine the minimum of the two minima next to the intensity maximum representing 
our blob, and subtract this value from /',' j . Lastly, we perform the gaussian fit of the resulting 
intensity profile over the interval /,. Several conditions must be met by the intensity profile 
and the resulting cr-error from the fit in order for it to be considered in the length calculation. 
An important check is that the variation in the intensity of the background profile I'^ ^^ 
along the interval /, is different than that of the blob's intensity profile. For that, we perform 
a gaussian fit of /^'^^ over /, and ensure that the resulting FWHM of the latter is smaller than 
the FWHM obtained from the fit of the blob's intensity profile, otherwise the fit is rejected. 
If the fit satisfies the test we proceed to define the length Z''^''''' of the blob as the FWHM 
resulting from the fit. The length /'•* of the blob at the i'-th measurement is then defined as 
the mean over the lengths Z'-* for each {x,t) e (x,',r,)i. A standard deviation for each i-th 
measurement of length is therefore also obtained, offering an estimate of the error involved. 
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The method for calculating the widths w''* is similar, with the difference that the gaussian fits 
are performed over the orthogonal direction to the path of the blob (right panel of Figure|2}- 



3.3 Temperatures 

The determination of the temperature T''^ for the blob is based entirely on the determination 
of its spectral profile as described in Section ISTI Given the blob's spectral profile /^'j, we 
assume that its broadening is mostly due to thermal motions. This assumption is supported 
by the fact that the observed plasma moves mostly one-dimensionally, (in the direction of 
the magnetic field, since it is partially ionized), and unresolved plasma motions are confined 
to the very small observed widths (transversal component to the magnetic field). However, 
turbulence and the Stark effect can be important, and thus the temperatures calculated here 
are only an upper limit. 

We thus assume that the FWHM'* obtained from the gaussian fit to /',' j is related to 
temperature according to: 



Aq llk^T 



FWHM'-'' =2V21n2 — ^ — S + v^ , (2) 

where Xq denotes the Ho- line center, inn is the hydrogen mass, k^ is Boltzmann's constant, 
c is the speed of light and v^mc is the microscopic velocity accounting, for instance, for 
turbulence. Assuming that this last term is 0, we obtain upper bounds for the plasma tem- 
perature. Since we allow the wavelength interval over which the gaussian fit is made to vary 
(as described in Section lSTt . we obtain an estimate of the errors involved in the temperature 
calculation. 



4 Results 

In this Section we present the results of the statistical analysis of the observed on-disk coro- 
nal rain blobs for the various datasets. The results presented here are obtained with the 
methods explained in Section [3] In total, 309 single coronal rain blobs were successfully 
traced resulting in 471 projected velocity measurements from the space-time diagrams ex- 
tracted along the paths drawn in the five data sets. By far the most coronal rain occurrences 
were measured in data set 1 (235 blobs), while in data sets 2, 3 and 5 on the order of 20-30 
blobs were found. Data set 4 only resulted in 2 measured coronal rain blobs, both following 
the same path, as can be seen from panels d in Figure [T] 



4.1 Dynamics 

Following the procedures outlined in Section lSTI we determined the slopes of the dark traces 
in the space-time diagrams in order to derive the projected velocities of the coronal rain 
blobs. These measurements were subsequently used as a starting point to determine Doppler 
component velocities, as well as the magnitude of the total velocity vector. Figure |2] shows 
the obtained velocity distributions. 

In the top left panel of that figure the total velocity distributions are shown, both for 
the cumulative sample (black solid line) and all five data sets separately (differentiated by 
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different color and line styles, see the caption for more details). In the top right panel the 
Doppler velocities are presented, where positive and negative values correspond, respec- 
tively, to blueshifts and redshifts. As a Doppler velocity is required in order to derive the 
total velocity of a blob and not all measurements result in a successful determination of the 
Doppler component velocity, the number of measurements for which a total velocity could 
be determined is smaller than the number of projected velocities and hence the cumulative 
sample is the same in size as in the top right panel (showing the Doppler velocity distri- 
butions), but smaller than in the bottom right panel. This is due to the narrowness of the 
spectral windows in these data sets. When this is the case, only part of the full spectral pro- 
file of the blob is observed. Consequently, the difference between the background spectral 
profile and the (partially) observed blob profile will give a profile whose maximum can very 
well not be in the considered wavelength range. The estimation of the Doppler velocity by 
gaussian fitting or first moment calculation (see Section 13.1b delivers too large errors and 
is therefore not included in the statistics. Therefore, the absence of high Doppler velocities 
in the top right panel of Figure |3] is not due mainly to the lack of them, but to the increase 
in errors in the calculation for those cases. Also, due to the relatively low cadence and the 
small field of view, high velocity blobs can only be observed over a couple of time steps, 
and are thus much harder to detect. Nonetheless, comparison between the top left panel and 
the bottom right panel shows that the overall distribution is shifted to higher velocities (from 
3-111 km s"' to 28-120 km s"'), as well as being slightly more skewed displaying a longer 
tail to higher velocities than to lower velocities. The average total velocity is shifted cor- 
respondingly upwards from an average projected velocity of 50 km s"' to 60km s"', with 
the averages for the several total velocity subsamples deviating typically by up to 20 km s"' 
from that value. However, it should be noted that the cumulative distribution peaks at a lower 
value than the average, between 45-55 km s"'. The lower velocities determined (below the 
average) are almost solely contributed by data set 1, as is the case for the projected velocities. 

The bottom left panel of Figure |3] shows the distribution of the total accelerations (i.e., 
accelerations derived from multiple increasing or decreasing total velocity measurements of 
one blob). Measurements from only three data sets (1, 4 and 5) were included in this figure 
as no accelerations could be properly measured in the other two data sets. The cumulative 
distribution displays a slightly longer tail to higher (positive) accelerations than to the de- 
celeration (i.e., negative acceleration) part and contains measurements with values in the 
range between -217 m s"- and 529 m s"^. The average value over the 156 accelerations is 
found to be 122 m s"-, while considering the accelerations (144 measurements) and deceler- 
ations (12 measurements) separately the averages are 137 m s"^ and -46 m s"^, respectively. 
Note however that also for the accelerations the distribution peaks below the average value, 
between 50-100 m s"^. 



4.2 Sizes 

Figure |4] shows the length and width distributions for those coronal rain blobs for which 
the size could successfully be determined according to the procedures described in Sec- 
tion [T2l resulting in 299 length and 309 width determinations. Overall, the lengths range 
between 217 km and 2.37 Mm with an average of 0.89 Mm, however the distribution peaks 
at a smaller size of around 600 km and the higher average is thus mainly a result of the 
extended tail above about 1.5 Mm. The histograms for the diflFerent data sets show a simi- 
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Figure 3 Velocity and acceleration distributions of the coronal rain blobs measured. Clockwise from the 
top left the panels show the distributions for the total velocity, Doppler velocity, projected velocity and the 
acceleration. In all panels the distribution for each data set has been diiferentiated by color and line style: 
the cumulative sample (black solid hne), data set 1 (red dashed hne), data set 2 (blue dotted line), data set 
3 (orange dash-dotted line), data set 4 (green dash-dot-dotted line) and data set 5 (purple long dashed line). 
The vertical dashed hne indicates the average value for the cumulative sample. The grey shaded bins in the 
top right panel indicate measurements for which the standard deviation exceeds 2 km s" ' . Note that in the 
lower left panel only accelerations are displayed for data sets 1, 4 and 5, and the coiresponding cumulative 
sample, as no accelerations could be determined for the other two data sets. The bin size for the upper panels 
and the bottom right panel is 5 km s"', while that for the bottom left panel is 0.05 km s"-. The total number 
of measurements £ N is indicated in the top right comer of each panel. 



lar behavior, having typical averages between 0.5-1.0 Mm and mostly extending only up to 
1.5 Mm. 

The cumulative width distribution, shown in the right panel of the same figure, ranges 
from 129-515 km, peaking around the same location as where the average of 324 km lies. 
Again the diff^erences between coronal rain blobs measured in the five data sets is small, 
which is reflected by the averages of the subsamples falling within 50 km around the cumu- 
lative average. 



4.3 Temperatures 



As described in Section 1331 the determination of the Doppler velocity by means of a gaus- 
sian fit also allows for the derivation of the blob's temperature. In accordance with Equa- 
tion ^ the profile width was thus converted to a temperature. Figure |5] shows the resulting 
distributions in a similar format as before, but now also including a grey shaded histogram 
corresponding to those measurements where the 1-cr value of the determined profile width 
divided by that same width exceeds 0.1. As indicated by this shaded histogram the trastwor- 
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Figure 4 Size distributions showing both lengths {left panel) and widths (right panel) for the traced coronal 
rain blobs. The bin size for the left panel is 200 km, while that for the right panel is 50 km. Format as for 
Figure |3] 



thiness of the temperatures in the lower bins is small, while most of the temperatures ob- 
tained from data sets 2, 3 and 5 actually fall in these lower temperature bins below 5x10^ K, 
where the cumulative distribution peaks. The average is found to lie at a higher value of 
7.8x10^^ K, which is a result of the high temperature tail reaching up to 22x10^ K. 
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Figure 5 Temperature distributions for the five data sets considered, as well as for the cumulative sample, 
with a bin size of 1.5x10' K. The grey shaded bins correspond to measurements where the l-o" value of the 
determined profile width exceeds 0.1 times that profile width. Further format as for Figure |3] 



5 Discussion 



Table [2] offers a comparison of the a verage values for the dynamics, sizes and temperatures 
between the coronal rain observed in lAntolin and Rouppe van der Voortl(l2012h (mostly off- 
limb) and the on-disk blobs observed here. Despite the large number difference in detected 
blobs, the resulting statistical distributions are very similar, thus suggesting the same phe- 
nomen on in both cases. Similar dynamics in coronal rain are also reported w i th other instru- 
ments iSchriivei 2001; De Groof. Berahmans. and van Driel-Gesztelvi 2004; De Groof et al.l 
12005; Miiller et al. 2005; Antolin. Shibata, and Vissers 2010; Antolin and Verwichte 20.U^ 
However, elongated blob-like cool and dense structures as those presented by coronal 
rain are not unique in the solar atmosphere and caution must be taken. Structures with simi- 
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Table 2 Average values for measured chai'acteristics between 
lAntolin and RoupDe van der Voorll (120121) ('off-limb' case) and the 
disk' case) 


the 
on-disk 


coronal rain reported in 
blobs observed here ('on- 


Case 


Total blob number Total velocities 
kms-' 


Accelerations 
kms-2 


Lengths Widths 
[Mm] [Mm] 


Temperatures 

[lO'K] 


Off-limb 
On-disk 


2552 67.1 
309 60.7 


0.0835 
0.1224 


0.74 
0.88 


0.31 
0.32 


9.6 

7.8 



lar features are being reported in prominences, as matter detaching from the main structure. 
Such observations take place mostly for quiescent rather than active prominences (or fila- 
ments], with hig h resolution instruments such as SST/CRISP, Hinode/SOT and SDO/AIA 
JChaee t al. 2Q0 8l:lLinl20l"ol:IZapi6r and Rudawvl2010l : lBerger et al.l2010l : lHm"ier, Isobe. and Watanabd 
I2OIIL iLiu. Berger. and Lowll2012h . Apart from eruptive prominences or cases of ejected 



blobs (where velocities are usually high), the reported velocities for these prominence blobs 
are often on the order of 30 km s"' , values that are usually much lower than the average val- 
ues for coronal rain blobs (in average 60-70 km s"'). A family of fas t moving prominence 



blobs with speeds between 40 and 115 km s ' has also been reported by van Noort and Rouppe van der Voord 
( l200a ) and more recently by Lin. Engvold. and Rouppe van der VoortI ( 120 12h in active 



re- 



gion prominences. These show constant traveling speeds along rather straight-like magnetic 
bundles with non-thermal velocities below 9 km s"'. Coronal rain can exhibit broad ve- 
locity distributions, including slow falling blobs with velocities as low* as 10 km s~' or so 
JAntolin. Shibata. and Vissersll20ld : lAntolin and Rouppe van der Voortll2012h . thus indicat- 
ing that a distinction criterion solely based on speed is incorrect. Since it is easy to discern 
prominences during observations, it provides a basis for distinguishing both kinds of blobs. 

When observing on-disk, structure characterization usually becomes more complicated. 
For our case, this is so since the condensations forming coronal rain are chromospheric 
in nature and the projection effect eliminates all obvious ways of calculating the height. 
Therefore, deducing the 'coronal' character of observed on-disk blobs is far from a trivial 
matter. The observed fast velocities in coronal rain can easily be explained due to the coronal 
heights at which the blobs form. Since the observed paths followed by the blobs are always 
loop-like (order of magnitude estimates indicate that the Ho- radiation coming from neutral 
Hydrogen in the rain is strongly coupled to the i ons component, and therefore that corona l 
rain is a good tracer of the coronal magnetic field. lAntolin and Rouppe van der Voortll2012h . 
the blobs have long paths over which the effective gravity along the loop can act. It is an 
established fact, however, that coronal rain blobs fall at velocities smaller than free fall. The 
calculated accelerations are accordingly lower than the effective gravity along loops, as is 
also observed here (bottom right panel in Figure |3). Therefore, it would seem that a fast 
falling cool and dense blob is most of the time a signature of coronal rain. Accepting it 
(generally) as a sufficient condition, it is however not a necessary condition. 

What really defines coronal rain is the sudden appearance in the corona (on the timescale 
of minutes) and subsequent fall along loop-like structures. The sudden appearance is char- 
acteristic of the thermal non-equilibrium mechanism. Also, coronal rain interpretation can 
be given if blobs are observed accelerating towards both footpoints of a loop (bidirectional 
motion), although the clear determination of the two legs of a loop can be extremely difficult. 
Further support can be found if the sudden appearance is accompanied by velocities that are 
at first very small (both projected and Doppler velocities), and subsequently increase grad- 
ually. This corresponds to blobs forming near the apex of loops, where the effective gravity 
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is small. Caution must be taken, however, since the projection effect and the line of sight 
angle with respect to the path geometry can also make a siphon flow (matter moving from 
one footpoint to the other) suddenly appear and gradually accelerate. Such a scenario can 
be rejected if the total velocity (calculated from the Doppler and projected components) is 
known along the trajectory. 

Apart from filament blobs, other on-disk blob-like features have been recently observed 
bv lVissers and Rouppe van der VoortI (120121) . Termed 'flocculent flows' by the authors, the 
blobs are observed close to a sunspot along short Ho- filaments with lengths around 15 Mm. 
The observed velocities show a mean around 35 km s"'with broad tails going as high as 
80 km s~', and seem to occur in a non-continuous way. In some cases, seemingly contin- 
uous motion of the blobs from one footpoint to the other is implied in the observations, 
suggesting that continuous siphon flow-like motions may be present. However, due to the 
blob-like shapes, the non-continuous occurrence character, the high velocity component and 
some bidirectional flows found, the authors argue that siphon flows cannot solely explain the 
phenomenon. Thus, a combination of both thermal non-equilibrium and pressure diff^erence 
between footpoints is suggested as a plausible explanation for flocculent flows. 

All together, the presented characteristics for our on-disk blobs differ from other seem- 
ingly non coronal rain-related blob-like structures observed so far, but do match closely, in 
a statistical sense, the coronal rain characteristics. We thus feel confident in concluding that 
the observed blobs correspond to the on-disk counterpart of coronal rain. It is interesting 
to note, however, that despite the high number of datasets analyzed, only 309 blobs were 
detected, which is a bare 12 % of the number of blobs detected oft-limb for one dataset (al- 
though a relatively long-lasting one). Although this fact could be taken as proof for a low 
frequency occurrence of the phenomenon, we believe that this is mainly due to the following 
reasons. First, the very small sizes involved are close to the diffraction limit of the telescope. 
Second, the not very wide spectral windows of the present datasets limit considerably the 
detection of high speed Doppler components for the blobs, which disables the calculation 
of the total velocities. Third, the cadence around 20 s in most of our datasets is significantly 
long for the short timescales involved in coronal rain, only allowing the detection of fast 
moving blobs on the small FOV available over a few time steps. And especially fourth, the 
fact that most blobs are very faint, and thus not very opaque when observed on-disk, leading 
to a very little intensity contrast with respect to the background intensity. These points make 
their on-disk detection extremely difficult. On the other hand, in the presented datasets we 
have two active regions at different times along their lives. The fact that coronal rain can be 
spotted at different times in an active region may be an indication of the common character 
of this phenomenon. 



6 Conclusions 

In this work several datasets corresponding to two active regions observed on-disk on differ- 
ent days (at different ju values) were analyzed in Ha with the CRISP instrument at the SST. 
On-disk blobs with special characteristics in all five datasets were selected and a statistical 
study of dynamics (total velocities, together with its Doppler and projected components), 
sizes (lengths and widths), and tem peratures was presented. Through compari son with the 
large statistical dataset presented in lAntolin and Rouppe van der Voorll (120 12h we showed 
that for all characteristics similar distributions exist. Other possible cool and dense blob-like 
features reported elsewhere were considered for the interpretation, such as blobs observed 
in prominences/filaments, or, especially, flocculent flows observed in the chromospheric 
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canopy bv lvissers and Rouppe van der Voord (l2012h . But neither match well the presented 
distributions shown here. The little number of blobs fo und, with respect to the large number 
presented in lAntolin and Rouppe van der Voord (120121) was also discussed, and considered 
to be most likely linked to the small and faint character of the phenomenon, which offers 
very little intensity contrast with the background. We have thus concluded that these on-disk 
blobs most likely correspond to the on-disk counterpart of the coronal rain, which was so 
far only observed off-limb. The fact that the active regions display coronal rain at different 
times during their life may be an indication of the common character of the phenomenon. 
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